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Abstract. We study the effect of an external field on (1 + 1) and (2 + 1) dimensional elastic manifolds, at 
zero temperature and with random bond disorder. Due to the glassy energy landscape the configuration of 
a manifold changes often in abrupt, "first order" -type of large jumps when the field is applied. First the 
scaling behavior of the energy gap between the global energy minimum and the next lowest minimum of the 
manifold is considered, by employing exact ground state calculations and an extreme statistics argument. 
The scaling has a logarithmic prefactor originating from the number of the minima in the landscape, 
and reads AEi ~ L* [ln(LzL~^)]~^''^, where ^ is the roughness exponent and 6 is the energy fluctuation 
exponent of the manifold, L is the linear size of the manifold, and Lz is the system height. The gap 
scaling is extended to the case of a finite external field and yields for the susceptibility of the manifolds 
Xtot ^ L'^^'^^~^ [{1 — In(L)]^'^^. We also present a mean field argument for the finite size scaling of the 
first jump field, hi ~ L"^"*. The implications to wetting in random systems, to finite-temperature behavior 
and the relation to Kardar-Parisi-Zhang non-equilibrium surface growth are discussed. 

PACS. 75.50.Lk Spin glasses and other random magnets - 05.70.Np Interface and surface thermodynamics 
- 68. 08. Be Wetting - 74.60.G Flux pinning, flux creep, and flux-line lattice dynamics 



1 Introduction 

In this paper we study zero-temperature or ground state 
elastic manifolds that are roughened by bulk disorder, in 
the presence of an external field. Such objects are rele- 
vant in many contexts of condensed matter and statistical 
physics 10,^. The essential point here is a competition 
between the elasticity, which prefers flat manifolds, and 
the disorder, which induces wandering in order to take 
advantage of the low energy regions in the system. This 
leads to a glassy, complicated energy landscape and the 
fact that the quenched randomness dominates thermal ef- 
fects at low temperatures. In two embedding dimensions 
(2D) such manifolds are under the name directed poly- 
mers (DP) [^^,|^ particularly interesting through their 
connection to the celebrated Kardar-Parisi-Zhang (KPZ) 
equation of nonlinear surface growth. Directed polymers 
have an experimental realization as vortex-lines in granu- 
lar superconductors Q . In higher dimensions elastic man- 
ifolds may be best considered as domain walls (DW) in 
ferromagnets with quenched impurities. Elastic manifolds 
have also other connections to charged density waves, the 
sine-Gordon model with disorder, random substrate prob- 
lems, and vortex lattices, to name but a few 

Let us start by introducing the classical spin-half ran- 
dom Ising Hamiltonian 



H 



I J Si Sj 



(1) 



where Jy is the coupling constant between the nearest- 
neighbor spins Si and Sj , and Hi is a field assigned to each 
spin. To this system we apply antiperiodic or domain wall 
-enforcing boundary conditions in one direction. The spins 
in the opposite boundaries, let us define in z-direction, z = 
and z = Lz are forced to be up and down, respectively. In 
the case of ferromagnetic random bond (RB) Ising systems 
one has Jy > and Hi = 0. In the minimum energy state 
the spins prefer to align on each side of the induced domain 
wall. When ^ the spins become frustrated and the 
task to find the ground state (GS) structure is most often 
related to spin glass physics On the other hand, 

when for simplicity = const. = J > 0, and Hi ^0 one 
arrives at random field (RF) Ising systems. The random 
field Ising model (RFIM) has an experimental realization 
as a diluted antiferromagnet in a field. In the RFIM the 
ferromagnetic couplings compete with the random field 
contribution, which prefers in the ground state to have the 
spins to be oriented towards the field assigned to them. 

Here we concentrate mainly on the random bond Ising 
Hamiltonian, i.e., Jy > and H = 0, and in some special 
cases extend the discussion to RF domain walls as well. In 
the simplest case the spins are located in a square lattice, 
in d = 2, or a cube in d = 3, so that the lattice orien- 
tation is in the {10} and {100} directions, respectively. 
The elastic manifold is the interface, with the dimension 
D = d—1, which divides the system in two parts of up and 
down spins. At T = the problem of finding the ground 
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state domain wall, which minimizes the path consisting 
of unsatisfied bonds between the spins on opposite sides 
of the domain wall becomes "global optimization" . In our 
case the displacement field is one-dimensional, n — I, but 
one can certainly think about generalizations so that the 
total dimension of a system d = {D + n), where n > I is 
the dimension in which the manifold is able to fluctuate. 

The continuum Hamiltonian of elastic manifolds with 
an external field and n = 1 may be written as 



r 



{Vz(x)}2 + F.{x, z(x)} + /i{z(x)} 



(2) 

The elastic energy is proportional to the area of the inter- 
face given by the first term, and F is the surface stiffness 
of the interface. The second term of the integrand comes 
from the random potential, and the last term accounts 
for the potential caused by the external field. The use of 
random bond disorder means that the random potential is 
delta-point correlated, i.e., (l/r(x, z')) = 2X'5(x — 

x')5(z — z'). In the random field case z)Vr{'x' ^ z')) ~ 

^(x — x')(z — z'). The Hamiltonian (||) is also applicable to 
wetting in a three phase system, where two of the phases 
are separated by an interface in a random bulk |2[|l|,|l|. 
In that case h{z) is equivalent to the chemical potential, 
which tries to bind the interface to the wall, and competes 
with the random potential, in the presence of which the 
interface tends to wander in the low energy regions of the 
system. 

Below the upper critical dimension the geometric be- 
havior of elastic manifolds is characterized by the spatial 
fluctuations. For the mean-square fluctuations one has 



z(x) - z(x) 



(3) 



where z(x) is the height of the interface with the mean z, 
x is the D — d — n dimensional internal coordinate of the 
manifold, L is the linear size of the system, and ^ is the 
corresponding roughness exponent. At low temperatures 
in {D + n) — {1 + I) dimensions with RB disorder, i.e., 
when one actually considers a directed polymer in random 
media the roughness exponent is calculated exactly 

via the KPZ formalism to be C = 2/3. In higher manifold 
dimensions D with n = I the functional renormalization 
group (FRG) calculations give the approximatively values 
C ^ 0.208(4 - D) for RB disorder and ( ^ {4- D)/3 
for RF disorder |j,|T|]. The expression for C. tells also that 
the upper critical dimension for the elastic manifold is 
I?u = 4. For manifolds with varying n and D Balents and 
D. Fisher have derived using FRG C — [(4 — D)/{n + 
4)]{l + (l/4e)2-[("+2)/2l[(n + 2)V(n + 4)][l-...]} |l§. At 
zero temperature the total average energy E of an elastic 
manifold equals its free energy and grows linearly with 
the system size and its fluctuations scale for all n as 

AE = {{E-Eff^ - L«, where 



2C + - 2 



(4) 



is the first non-analytic correction to the energy. The same 
hyperscaling law holds for RF manifolds, too, in D > \ 



p8[ . Having a positive 6 implies that the temperature is 
an irrelevant variable in the renormalization group (RG) 
sense and the T — Q fixed point dominates. For 13 = 1, 
n > 2 there exists a Tc, and T = fixed point domi- 
nates only below Tc, likewise always for n < 2 ||T^. At 
the randomness dominated pinned phase the temperature 
is "dangerously" irrelevant, which means that in RG cal- 
culations the interesting correlation functions cannot be 
obtained by setting T to zero. Above Tc the fluctuations 
become random walk -like with C, = 1/2 and — Q. For 
D — 1 with n > 1 there is no exact result existing for the 
roughness exponent below and hence whether — > 1/2 
and — > for a finite ric is still an open question - that is, 
what is the upper critical dimension of the KPZ equation. 

Elastic manifolds self-average in the sense that the in- 
tensive fluctuations of the roughness and the energy de- 
crease with system size. However, in this paper we will 
show that, due to the glassiness of the energy landscape, 
for example the behavior of the mean position of a typ- 
ical example of a manifold z(x) does not coincide with 
the disorder average, here denoted by (), over many real- 
izations with different random configurations. Introducing 
the external field to a random system induces often drastic 
changes and one experimental possibility is to study the 
ground state behavior e.g. by measuring the susceptibility. 
In disordered superconductors the external field is due to 
the current density j, which drives the vortex lines. 

Here we study the elastic manifolds in an external field 
while being especially interested in the information it gives 
about the energy landscape of the random system [ p0| , 
which is closely related to the susceptibility. The external 
field is appHed to the coupHng constants [see Eq. (||)] so 
that J±(z) = J random + h(z) , whcrc J_L are ones in the 
z-direction and h is the amplitude of h{z), the strength of 
the external field. Since we have the field potential linear 
in height, h{z) = hz, in ferromagnetic random systems 
with a domain wall the external field h{z) may be trans- 
formed to a constant external field term —H Si in the 
Hamiltonian (0) . Our results are a systematic extension of 
early numerical work by Mezard and relate to more recent 
discussion of the energy landscapes of directed polymers 
by Hwa and D. Fisher |gl|,^. Due to the glassiness of the 
energy landscape the position changes in "first order type" 
large jumps, at sample-dependent values of the external 
field. The second scenario in which the perturbations take 
place locally via "droplet" -like excitations is ruled out by 
a scaling argument and numerical results. 

We start by studying a specific case where the number 
of valleys is fixed to a constant. We first consider the case 
without the external field, and find the scaling behavior of 
the energy difference of the global energy minimum and 
the next lowest minimum of a manifold numerically and 
derive it also from an extreme statistics argument. The 
scaling of the gap does not only depend on the energy 
fluctuation scale, but has in addition a logarithmic fac- 
tor dependent on the number of the low energy minima. 
The gap scaling is extended to the finite external field 
case, to derive the susceptibility of the manifolds under 
the assumption that the zero-field energy landscape is still 
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relevant. The corresponding numerics allows us also to de- 
duce the effective gap probability distribution without any 
a priori assumptions. We then proceed by constructing a 
mean field argument for the finite size scaling of the first 
jump field in a more general case. This agrees with "grand- 
canonical" numerics in which the manifold is allowed to 
have an arbitrary ground state position in the system: 
the number of valleys in the energy landscape fluctuates 
and the important physics is caught by the simple scal- 
ing considerations. Finally, we discuss wetting in random 
systems in the light of our results, and compare with the 
necklace theory of M. Fisher and Lipowsky that applies in 
the limit of large fields. We also consider finite tempera- 
tures, and relate the physics at finite external fields to the 
Kardar-Parisi-Zhang growth problem through the arrival 
time mapping. The implications of our results in that case 
concern the two first arrival times and their difference of 
an interface to a prefixed height. Note that we have pub- 
lished short accounts of some of the work contained here 
earher|2|,|||. 

The paper is organized so that it starts with a short 
review of the thermodynamic behavior of the interfaces in 
random media at zero temperature in the presence of an 
applied field, and of the folklore concerning energy land- 
scapes, in Section ^. Section ^ introduces the mean-field 
level behavior of the interface when the external field is 
applied. The analytical probability distribution of the first 
jump field is derived and its relation to the susceptibility 
is discussed in Section ||. In the section also the lowest en- 
ergy gaps between the local energy minima and the global 
minimum are derived from extreme statistics arguments. 
The numerical method of calculating the first jump field 
hi of an interfac e in a fixed height is introduced in Sec- 
tion |]^. Section |5.2| contains the numerical results of the 
first jump field in (I + I) and (2 -|- 1) dimensions with the 
corresponding jump geometry statistics as well as the en- 
ergy gaps of the first lowest energy minima. The results 
are compared with the analytical arguments presented in 
the previous section. In Section ^ a mean field argument 
for the first jump field of an interface which lies originally 
in an arbitrary height is derived and compared with the 
numerical data and the evolution of the jump size distri- 
butions is studied. As an application for the finite field 
case the wetting phenomenon is considered and the cor- 
responding wetting exponents are numerically studied in 
Section ^. The arguments presented in the paper are dis- 
cussed from the viewpoint of int erfac e behavior at low 
but finite temperatures in Section 8.1 and in th e co ntext 
of first arrival times in KPZ growth in Section ^.2\ The 
paper is finished with conclusions in Section ^. 



2 Thermodynamic considerations 

The interesting behavior arises since in the low tempera- 
ture phase e.g. directed polymers are found to have anoma- 
lous fluctuations resulting from the regions of the random 
potential with almost degenerate energy minima, which 
are separated by large energy barriers | 



9,21 



, These 

spatially large-scale low-energy excitations are rare, but 



are expected to dominate the thermodynamic properties 
and cause large variations in the structural properties at 
low temperatures. We will use the arguments in the next 
section in order to study the movement of the manifolds 
in equilibrium when an external field is applied. Hwa and 
D. Fisher derived that for a polymer fixed at one end 
a small energy excitation with a transverse scale A 
where / is the linear length of the excitation), scales as 
Z\^/^ and gives rise to large sample-to-sample variations in 
the two-point correlation function and dominates the dis- 
order averages. The probability distribution for the sizes 
of the excitations was found to be a power-law from the 
normalization of "density of states" of small energy exci- 
tations W - Z\-"r^ - j^/^«+2-i/c^ gin(,g 6* = 2C - 1, 
Eq. ^ . This is in contrast to the high temperature phase 
one, where W is Gaussian. For the finite field case Hwa 
and D. Fisher argued that the power-law distribution for 
the excitations is due to the statistical tilt symmetry, 
which means that the random part of the new poten- 
tial of the polymer is statistically the same as the old 
one, when the polymer is excited by an applied field. In 
Ref. [0 Mezard showed numerically that the energy dif- 
ference of two copies of the polymer in the same realiza- 
tion of the random potential scales as , where 6 = 1/3 
in d = 2. He also showed that the ground state configura- 
tion of a polymer, which is fixed at one end and applied 
with an external field h to the end point of the polymer, 
i.e., h{z{x.)} = hz5{x — L) in Hamiltonian (||), changes 
abruptly with a distance i^, C = 2/3, when the field is in- 
creased by . Calculating the variance of the end point of 
the polymer var{z) ^ (z^) — (2)^ Mezard found that typ- 
ically it is zero, i.e., the ground state configuration does 
not change. With the probability L^^ , there is a sample 
for which /i = is the critical value of /i, and the config- 
uration changes with a large difference compared to the 
original state L^^ ^ Eq. (Q). Hence the total av- 

erage variance becomes var{z)ave ~ L. In this paper we 
generalize such studies as discussed later. 

The physics of DP's have been shown by Parisi to obey 
weakly broken replica symmetry , a "baby-spin glass" 
phenomenon. This means that the replica symmetric so- 
lution of a DP is degenerate with the solution with the 
broken replica symmetry. This is also evidence for an en- 
ergy landscape with several, far away from each other, 
nearly degenerate local minima. Mezard and Bouchaud 
studied later the connection between extreme statistics 
and one step replica symmetry breaking in the random en- 
ergy model (REM) p6| , in a finite, one-dimensional form 
[ P7[ . They computed using extreme statistics the prob- 
ability distribution of the minimum of all the energies 
-\- E{x), where E{x) are random and follow a suit- 
able probability distribution so that its tails decay faster 
than any power-law, e.g. Gaussian. This can be seen as a 
toy model for an interface in a random media, in which 
case the quadratic part plays the role of the elasticity. 
The minimum of the total energy is easily seen to be dis- 
tributed according to the Gumbel distribution |28 of ex- 
treme statistics, and likewise for the position of the in- 
terface the shape of the probability distribution to be ap- 
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proximately Gaussian, fn Section ^ we will use extreme 
statistics to study the scaling of the lowest minima and 
energy gaps in elastic manifolds in a situation analogous 
to that of Bouchaud and Mezard. 

The response of the manifolds pinned by quenched im- 
purities to perturbations was discussed in terms of ther- 
modynamic functions by Shapir pof , including the sus- 
ceptibility of the manifolds. The formal definition of the 
susceptibility for a D dimensional manifold in a d dimen- 
sional embedding space reads: 



/ dm\ 
h^o+ \ oh / 



(5) 



where the change in the magnetization of the whole d di- 
mensional system is calculated in the limit of the vanishing 
external field from the positi ve si de and the brackets imply 
disorder average. In Section 5.2 we derive the susceptibil- 
ity based on energy landscape arguments. By assuming 
smooth, analytic behavior in the manifolds' thermody- 
namic functions, Shapir found that the susceptibility (for 
a surface of dimension D) is proportional to the displace- 
ment of the manifold in the limit of small field, which 
is applied uniformly to the whole manifold, and xd ~ 
(PE/dK^ ~ L''+2". a results from an argument concern- 
ing the energy gap and the external field: AE ~ hL'^L^ , 
where the external field couples to a droplet of size and 
the area of the manifold is . This should be equal to the 
energy difference: AE ~ L^+", and thus a = 2 — C,, since 
6 = 2( + D — 2. Hence the susceptibility was derived to be 
Xd ^ L^"*"^ and the susceptibility per unit hyper-surface 
vary as for a manifold of scale L. This surprisingly is 
independent of the type of the pinning randomness. 

The landscapes of DP's have also been studied by Jogi 
and Sornette Q by moving both end points step by step 
(with periodic boundary conditions) from (0, z) and (L, z) 
to (0, z+1) and {L, z + 1), and finding the ground state at 
each step. They studied the sizes of the avalanches as in 
self-organized systems, i.e., the areas the polymer covers 
when moving to the next position. They found a power- 
law for the avalanche sizes P{S) dS ^ dS^ where 
m = 2/5, up to the maximum size Smax ~ L^^^- Smax is 
actually the size of the largest fluctuation of the polymer 
in such a controlled movement, i.e., L x L^. The authors 
did not however study the energy variations of the mani- 
fold during the "dragging process" . 

In the limit of a large h the Hamiltonian of Eq. (||) 
is applicable to (complete) wetting interface in random 
systems with three phases. In this case the manifold is 
the interface between i) non-wet and ii) wet phases near 
iii) a hard wall at z = 0. This problem was studied by 
mean-field arguments by first Lipowsky and M. Fisher 
and later elaborated numerically The choice of 

wetting potential that corresponds to the linear field used 
here is called in the wetting literature the weak-fluctuation 
regime (WFL). In the WFL regime one can derive the wet- 
ting exponent -0 by a Flory argument, since the mean dis- 
tance of the interface from the inert wall z is of the order 
of the interface transverse fluctuations see Fig. |l|, i.e., 
the field is strong enough to bind the interface to the wall. 



The Hamiltonian can now be minimized by estimating the 
total potential to be Vtot = Vfi (z) + Vw ("z) • ~ hz is 
the wetting potential induced by the external field h. The 
fluctuation-induced potential Vfi is the potential of the 
free interface 14; + = ■j{V2:(x)}^ -|- K(x, z), which 
gives, when minimized, ~ "^^^ elastic energy term 
V^i gives the scaling of Vfi with respect to the correlation 
lengths, and expanding the square-term of the gradient of 

the height gives const i-|-const2 ^ -I- higher order terms. 

So, Vji ~ ^J"^, where the decay exponent t = (2 — 2^)/^ 
defines the scaling between the pinning and elastic ener- 
gies. Minimizing Vtot in WFL, where z ~ and is much 
larger than a lattice constant, gives 



-1/(1+^) 



Thus the wetting-exponent t/j = 1/(1 + t) becomes 



c 



2-C 



(6) 



(7) 



Huang et al. |^,|T3[ did numerical calculations for directed 
polymers at zero temperature using transfer-matrix tech- 
niques for slab geometries, L ^ L^, where L is the length 
of the polymer, and is the height of the system, con- 
firming roughly the expected exponent ^ ~ 0.5. Our re- 
sults, presented in Section Q are in line with the expec- 
tations in (1-1-1) and we also present the first studies in 
(2-1-1) dimensions. 



3 Glassy energy landscapes: 
perturbations 



length scale of 



We next employ a version of the scaling arguments by Hwa 
and Fisher, and Mezard for the next optimal position of 
the directed polymer. The goal is to investigate the pre- 
ferred length scale of the change that takes place, i.e., the 
size of the "droplet" created with the field. The numer- 
ical calculations presented in this paper are done such a 
way, that first the ground state is searched and the con- 
figuration stored. After that the external field is applied 
and the energy is minimized again. If there is a change in 
the configuration, the difference in the mean height of the 
manifold, the so called "jump size" is analyzed. Most of 
the studies are done in isotropic systems, i.e., the height 
Lz of the system is of the order of the linear length of the 
manifold, Lz oz L. One finds that the behavior of the poly- 
mer has a first order character, the optimal length scale is 
such that the whole configuration changes, and this holds 
for higher dimensional manifolds, too. 

Hwa and Fisher defined the next optimal position of 
the DP fixed at one end and with displacement A to scale 
a,s A ^ L'', where L is the length of the polymer. Now one 
assumes that the energy gaps between two copies of the 
polymers with overlap, and energies Eq and Ei, grows as 
El — Eq ^ , i.e., is just energy fluctuation exponent of a 
polymer (as demonstrated by Mezard numerically). Then 
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as a generalization with n = 1 it follows that £'i — £'o ~ 
Z\^/^. The external field has a contribution for the energy 
differences of interfaces with the dimension D Ei — Eq ^ 
hL^A ~ ft,Z\^+-°/'> . Assuming that this difference balances 
the gap, and using the hyperscaling law for the roughness 
and energy fluctuation exponents 9 — + D — 2, we get 

The exponent a is negative assuming that the roughness 
exponent is below two, which is of course satisfied in the 
case of both RF and RB domain walls. Hence, the smaller 
the field the larger the excitations and thus large excita- 
tions are the preferred ones, at least below the upper crit- 
ical dimension. Besides RB systems Eq. (^) works only for 
RF systems with a ferromagnetic bulk (among others, the 
relation Ei ~ Eq ^ has not been demonstrated to hold 
for RF interfaces). If the bulk of a RF magnet is param- 
agnetic, the stiffness of DW's is expected vanish on large 
enough length scales, translational invariance is broken, 
and AE = Ei — Eq does not scale. The scaling relation 

should hold also for general n. One should note in the 
case h is applied in one direction but the droplet may ex- 
tend in n dimensions, in order to Eq. (|^) to hold, also the 
projection of the droplet in the applied field direction has 
to have the scaling as A^^'^. 

In Fig. ^ it is shown what happens for two different 
random realizations of DPs, when a perturbing external 
field is applied. Fig. ^(a) shows the mean height of the 
polymer normalized by its original height 'z/zq, when the 
external field is increased. The first realization 1° shows 
a large jump of a size half of the height of its original po- 
sition at first jump field hi = 8 x 10~^. In Fig. |2|(b) the 
two positions of 1°, before and after the first jump, zo{x) 
and Zi{x), are shown. The other scenario would be that 
the directed polymer would continuously undergo small 
geometric adjustments and get meanwhile bound to the 
wall, z = 0. An example of such a droplet, is shown in 
Fig. ||(b) as the second realization 2°, see also Fig. |l[ How- 
ever, a further field increase, demonstrated in Fig. ||(a) 
for 2°, leads to a global jump too. This leads to a pic- 
ture in which assuming a starting position far enough 
from the system boundary, a finite number of large jumps 
exists from the original position 2;o(x) to the positions 
zi(x), 22 (x), . . . , Zri(x), . . ., closer and closer to the wall, 
i.e., zo(x) > zi(x) > Z2(x) > . . . > z„(x) > 0. 

The mean jump length is defined to be Azn — "Zn-i — 
z„, and the jumps take place at the fields /ii, /i2, . . . , ^n, . . .. 
The field value hi corresponds the jump from zo(x) to 
zi(x). Finally, after a finite number of jumps, with the 
interface being in the proximity of the wall, around h ~ 
2 X 10~^ in Fig. ^(a), see also Fig. |l|, the interface evolves 
inside the last valley next to the wall continuously, and 
the wetting behavior, Eqs. and (Q), applies. We dis- 
cuss this picture of consequent jumps before the wetting 
regime in Section ^ together with the numerical results. 

The global changes (large jumps) induce finite changes 
in the magnetization (m), and are reminiscent of first or- 
der phase transitions. Note that the field hi of the first 
change or jump obeys for any particular ensemble a prob- 



ability distribution, and thus a co-existence follows be- 
tween systems that have undergone a jump or change in 
m and those that are still in the original state. This "tran- 
sition" is a result of level-crossing between the valleys in 
the energy landscape for the interface, sketched in Fig. ||. 
To change between the geometrically separated minima, 
an external field is applied, which at jumps plays the role 
of a latent heat. Originally the interface lies at height 
and has an energy Eq. When the field is applied inter- 
face's energy increases by hL^'zQ and at hi the interface 
jumps to li < Jq having energy E(hi) = Ei -\- hiL^'zi = 
Eq + hiL^zo, where Ei is the energy of the interface 
at zi without the field. Similar behavior takes place at 
hn, n — 2,3, . . ., when the interface moves from z„_i to 
Zri < "Zn-i, until the cross-over to the wetting regime. 

In the thermodynamic (TD) limit, i.e., for very large 
systems, when L^L^ —> oo, the interface may originally be 
located anywhere in the system. Assuming that the rough- 
ness exponent ^ does not define only the width of the man- 
ifold, but also the width of the minimum energy valley, i.e., 

is the only relevant length scale in the transverse di- 
rection, gives that the manifold should have ^ L^/L'^ 
minima from which to choose the global minimum posi- 
tion. For C < 1 the number of the minima grows with 
system size, if the geometry is kept isomorphic, (x L. 
The requirement C < 1 holds for all RB interfaces; and 
for RF ones with ferromagnetic bulk, when D > 1. In the 
special case that the local minimum closest to the wall 
coincides with the global minimum, the wetting regime 
is entered at once without any jumps. The prediction of 
Eq. (||) that large scale excitations should dominate is 
an asymptotic one, and for finite system sizes one may 
find also small excitations, which are less costly in energy 
(Fig. H). However, the fraction of these cases decreases 
with system size, see Fig. ^. It depicts in d = 2 up to 

= 1000^, with at least N = 3000 realizations per data 
point, the probability of finding a non-zero overlap q be- 
tween the interfaces before and after the first jump, q is 
the disorder average of the fraction of the samples, which 
have overlap between the cases before and after the jump, 
i.e., for at least one x zq{x) = zi{x). The probability goes 
as g ~ 0.56L~°'^'^^*' °^, and confirms the expected behav- 
ior that for isotropic systems in the L — > oo limit the first 
jump is always without an overlap, i.e., a macroscopic one. 
A first guess would give that q ~ ~ L^'^ A similar 

trend seems to exist in (2 -I- 1) dimensions, too, but our 
statistics is not good enough to determine the exponent 
of q reliably. This is since q is averaged over a binary dis- 
tribution (overlapping and non-overlapping jumps), and 
thus very large ensembles would be needed. 

4 Extreme statistics arguments 

Next we compute the field hi and the susceptibility, Eq. (||), 
by taking into accont the "co-existence" of the original 
ground state and the state after the jump. Assuming that 
the relevant process in the response of a domain wall is 
the large-scale jump, the susceptibility per spin of a sys- 
tem with a domain wall, Eq. (|5|), may be written in the 
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form 



X 



h^a+ \ Ah / 



lim P(/»i), (9) 



where P{hi) is the probability distribution of the first 
jump fields with the corresponding first jump size Azi, 
and the magnetization m{h) ~ z(h)/Lz. The limit h 
0+ is taken over the probability of having a jump, and thus 
the susceptibility will reflect the co-existence phenomenon 
discussed above. One uses also a plausible assumption, 
found to be true in the numerics, that the change in the 
interface does not depend on the threshold field hi . A fur- 
ther simplification is obtained by assuming (this will be 
shown numerically later) that the jumps have an invariant 
size distribution, independent of L and L^, after a normal- 
ization with Lz- This gives (^^^i) ~ Lz- Hence the finite 
size scaling of the susceptibility per spin depends only on 
the probability distribution of the first jump field in the 
limit of the vanishing external field, and in particular on 
the "rare events" measured by such a distribution. The 
dependence of hi on the anomalous scaling is found to be 
true not only for directed polymers but also for the higher 
dimensional case. Such behavior is contrary to Shapir's 
work that assumed smooth, analytic thermodynamic func- 
tions and it is thus no surprise that th e sc aling of the 
susceptibility, derived in detail in Section 5.2, differs from 
that. 

Next we first derive the distribution for the first jump 
field and its relation to the lowest energy gap probability 
distributions. Then the probability distribution and the 
finite size scaling of the lowest gaps and also the lowest 
energy level are derived using extreme statistics arguments 
(see also [^2|,33 and our shorter account of parts of this 
work |2^). These gaps have a logarithmic dependence on 
the number of local minima. 

In order to derive the probability distribution of the 
first jump field P{hi), let us first calculate the probability 
Pq, that with a certain test field h' < hi, the interface has 
not changed yet, i.e. 



Po{hi > h') = 1 - Po{hi < h') = 1 - P 

Jo 



{h)dh. (10) 



By differentiating Po{hi > h'), one gets the probability 
distribution of the first jump field, since 

P{hi) = -^Po{hi<h')\„^,^. (11) 

We assume all the minima to be non-correlated and well 
separated from each other, see Fig. I|(a). There we have 
a global energy minimum Eo at zo (note, that here the 
energy values as well as the field contributions are nor- 
malized by , constant with a fixed system size), and 
energy gaps AEi — Ei—Eq and AEi* — Ei*—Eq > AEi 
with energies Ei at zi and Ei-^ at zi*, respectively. Then 
we apply the field h, i.e., we tilt the energy landscape, see 
Fig. ||(b). Due to the statistical tilt symmetry the external 
field h is assumed not to change the shape of the random 
landscape. At the lowest tilt to move the interface, i.e., at 



field hi, it moves to zi., since with the field hi £'0(^1) = 
Ea + hizo = Ei,{hi) = Ei^+hizi^ < Ei{hi) = Ei + hizi. 
However, sometimes Ei and Ei* coincide, so that the in- 
terface jumps to the true second lowest minimum. Now 
we index the positions of the Nz minima as Zi (need not 
to be in order), so that /i(zo ~ -^i) — hAzi, and we have 

Pa{hi > h') = ^^m\{AE, > hAz,) 
1=1 

= l[y-J^ mAE,)d{AE,)\. (12) 

By assuming all the gap energies from the same probabil- 
ity distribution P{AE), setting Azk^^ and taking the 
continuum-limit we get for the probability distribution of 
the first jump field using Eqs. (M), (111]) and (|l|). 



Pihi 



exp 



1 Jo 



P{AE) d{AE) dk 
P{khi/Nz) 



1- r 



P{AE)d{AE) 



dfc(13) 



There are two approaches to go on from the distribution 
of Eq. ([l3|): one may either compare its prediction with 
numerics with a trial distribution for P, or alternatively 
proceed by a complete Ansatz for the valley energies, from 
which the gap naturally follows. We now attempt the lat- 
ter, and comment on the first one together with the nu- 
merics. 

In order to calculate the probability distribution for 
the gap energies P{AE), we use a probability distribu- 
tion for the energy minima suitable for directed poly- 
mers, which has the advantage of decaying faster than any 
power-law, 'P{E) ~ exp{—BE^), ij > 0. For directed poly- 
mers with the "single valley" boundary condition case, i.e., 
one end of the manifold fixed and thus no room for other 
minima, it is known mostly numerically that the bulk of 
the distribution of the energy of a directed polymer V{E) 
is Gaussian (the exponent = 2 in the exponential), ex- 
cept for the tails, {E < Emm and E > Emax) in which 
different cut-offs, (?7_ = 1.6 and 77+ = 2.4) take over [|4[ 
1 1 . For the other manifolds that one might be interested in 
lie actual distributions have not generally been studied. 
In the following we usually approximate the distribution 
with Gaussian distribution (this is valid for 1/Nz > Emm 
and < Emax)- The validity of the approximation is dis- 
cussed in Section |5.2| , when the analytical arguments are 
compared with the numerics. 

The distribution [rj — 2 for Gaussian) for V{E) is writ- 
ten in the form: 



t 



V{E) = fcexp <^ - 



\E-{E)\ 
AE 



(14) 



where {E) ^ is the average energy of the manifold 
and k ~ {AE)^^ ^ normalizes the distribution. The 
extreme statistics argument goes so that in a system with 
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Lz/ L'^ minima the probability for the lowest energy to be E is 



LnAE)=N,V{E){1~Ci{E)} 



(15) 



where Ci{E) = V{e) de is called the error-function when rj — 2. For Gaussian V{E), likewise for general f] > 0, 



L]\[^(E) is known to be Gumbel distributed, exp(u — expw) 
distribution, Gn^{'^Ei, E) is given by 



The gap AEi follows similarly as Ln^E). Its 



GnA^Ei,E) = ^^^^ — ^-^V{E)V{E + AEi) 



{I - Ci{E + AE^)} 



(16) 



Gn^{AEi, E) is the probability that if the lowest energy manifold has an energy E, then the gap to the next lowest 
energy level is AEi. G^^iAEi, E) is actually a generalization of Gumbel distribution. Integrating Eq. ( p^ ) over all 
energies gives the probability distribution for the AEi 



PiAE,) = 

This becomes for Gaussian {rj = 2) and AEi <C (E) 



GNA^Ei,E)dE. 



(17) 



P{AE^) = 



N,{N,-l) 



k exp 



-{E - {E)f -2AEi{E - {E)) 
AE2 



k crfc 



E + AEi - (E) 
AE 



where erf denotes the error- function and erfc = 1— erf. Neglecting all O 



{fcerfc(^±^^)}^ 



terms and using a Taylor expansion around AEi — 0, i.e., Maclaurin-series, gives 



dE, (18) 
and higher order 



PiAEi) - 
1-{N,- 2)A:erf 



N.jN, - 1) 
2 

E-{E) 
AE 



k exp 



-{E-{E)y-2AEi{E-{E)) \ 



(N, - 2)AEikcxp 



AE^ 

-{E-{E)Y 
AE^ 



dE. 



(19) 
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The first two terms within the second parenthesis dom- 
inate for AEi small, and thus we see that to first order of 
P{AEi) ^ const (compare with the numerics presented 
below). In particular, one should notice that the probabil- 
ity distribution does not vanish for AEi ~ 0. 

This approach is very similar to the extreme statis- 
tics calculation of the one-dimensional version of the ran- 
dom energy model by Bouchaud and Mezard |^ (see 
Section ||). They derived the probability distribution for 
the location that gives minimum energy for the system, 
with the difference to our case that Hamiltonian reads 
+ E{x). The result becomes such that the distribu- 
tion for the position is approximately Gaussian, too. The 
one-dimensional system is close to our example except 
for the functional form of the external potential which is 
quadratic instead of a linear one. There is however one es- 
sential difference in that in their analysis the "field value" 
is fixed, whereas in our case we are interested in what 
happens in any particular sample as the field is varied. 
Nevertheless such a calculation could be compared to the 
distribution of the interface locations at a fixed field h. 

More important than the actual distributions is that 
the finite size scaling of the gap energies may be computed. 
Let us start by calculating the typical lowest energy level. 
The average of it is given by 



{Eo)= / ELNjE)dE, 



(20) 



and the typical value of the lowest energy is estimated 
from 

N.jniEo)) « 1. (21) 



Note, that in approximating the integral, Eq. (^0|), with 
the aid of the distribution, Eq. ( pT| ) (in the limit Ci [E) in 
Eq. (|r^ ) is small) the normalization 1/k should be taken 
into account. Eq. (|l|) gives 



(Eo) ^ {E)-AE[\n{N,)] 



l/r, 



(22) 



where AE ^ and for Gaussian t] = 2. 

To estimate the typical value of the gap, we make a 
similar approximation as in Eq. (|2l|) for {Eq), 



NAN. ~ 1) 
2fc2 



V{{Eo))V{{Eo) + {AE,)) « 1, 



(23) 



which with (Eg) and the fact that \ {AEi)\ ^ \{Eo)\ yields. 



{AEi) 



AE^f 



AE 



r,{{E) - {Eo))^-^ ,7 [ln(Ar,)] 



(24) 



We thus find that the gap scales as AE, - Z\£'[ln(7V^)]~i/2 
L^[ln(L2L^'')]^^/2, when rj — 2 and assuming as in the 
previous section, that ^ L./L''. If the interfaces are 
flat, i.e., C = which is true for {lOOlRB interfaces be- 
low the cross-over roughening scale Lc if ran- 
domness is weak, the same arguments should hold. How- 
ever, then the energy distribution is pure Gaussian, rj — 2, 
and 6 = D/2 due to Poissonian statistics. 



between the spins correspond to flow 
from a site St to its neighboring one 
The graph-theoretical optimization algorithm, a 



5 Susceptibility of manifolds 



5.1 Numerical method 



For the numerical calculations the RB Hamiltonian (||) 
with Jij > 0, H — 0, has been transformed to a random 
flow graph. The graph is formed by the Ising lattice and 
two extra sites: the source and the sink; and the coupling 
constants 2 J^^ 
capacities Cij 

I . 

maximum-flow minimum-cut algorithm, enables us to find 
the bottleneck, which restricts the amount of the flow that 
can get from the source to the sink given the capacities 
in such a random graph. This bottleneck, a path P which 
divides the system into two parts (sites connected to the 
sink and sites connected to the source) is the minimum 
cut of the graph and the sum of the capacities belonging 
to the cut Cij equals the maximum flow, the smallest 
of all cut-paths in the system. The source is connected to 
the sites in the Ising lattice which are forced to be up and 
the sink is connected to the sites which are forced to be 
down. The value of the maximum flow is the total min- 
imum energy of the domain wall, equivalent to the min- 
imum cut. The maximum flow algorithms can be proven 
to give the exact minimum cut for all the random graphs, 
in which the capacities are positive semidefinite and with 
a single source and sink The best known maximum 
flow method is by Ford and Fulkerson and called the aug- 
menting path method |Q. We have used a more sophis- 
ticated method called push-and-relabel by Goldberg and 
Tarjan which we have optimized for our purposes. It 
scales almost linearly, 0{n^''^), with the number of spins 
and gives the ground state DW in about minute for a mil- 
lion of spins in a workstation. Notice that one could use 
for DP's the usual transfer matrix method as well, but 
the max-flow implementation is convenient in the case of 
systematic perturbations to each sample. 

In this study we have done simulations for (1 + 1) 
and (2 + 1) dimensional manifolds. The system sizes ex- 
tend to L X = 1000^ and L'^ x = 400^ x 50. The 
number of realizations N ranges from 200 to 2 x 10"*. 
The random bonds are such that in 2D e [0, 1] uni- 
form distribution and Jij^x = 0-5 unless otherwise men- 
tioned. In three dimensions we have either 6 [0, 1] 
uniform distribution in all x,y, and z directions or dilu- 
tion type disorder, P{Jij) = pS{Jij — 1) + (1 — p)S{Jij), 
i.e., a bond has a value of unity or zero with the proba- 
bility p. We have used p — 0.5 and p = 0.95. The linear 
field contribution is applied in the z-directional bonds as 
J±{z) — Jrandorn + hz. When the jump field values with 
-'the corresponding jump distances are searched, the pre- 
cision is as small as Ah = 10~^ in order to be sure that 
no smaller changes would occur between the jumps, which 
could be the case if the precision was much larger. Periodic 
boundary conditions are used in x and y directions. 

When studying the susceptibility, we have controlled 
the number of the minima in the systems to avoid fluctu- 
ations in Nz (and the grand-canonical ensemble). To fix 
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the number of the minima {N^) ^ Lz/ in a certain sys- 
tem size, we have set the initial position of the interface 
zo in a fixed size window at height 'zq/ Lz — const. If the 
ground state interface is originally outside the window, 
large enough only for a single valley, it is discarded so 
that the statistics quoted are based on the successful at- 
tempts. If the window is well-separated in space from the 
2-directional boundaries it is obvious that this sampling 
has no effect on the statistical properties as the average 
energy. After an original ground state is found, with an 
energy Eq, the external field h is applied by increasing 
the cou plin gs by constant steps until the first jump is ob- 
served Q. We have also calculated the energy gap be- 
tween the global minimum and the next lowest minimum. 
In that case the initial position of the interface zq is also 
set to be in a fixed size window at height zq/L^ ~ const 
by discarding all the samples with the global minimum 
outside the window. After that the lattice is reduced so 
that bonds in and above the window are neglected and 
the new ground state, its energy E'l, and the correspond- 
ing gap energy AEi are found. Although the discarding 
procedure is slow we have at least N — 500 realizations 
up to system sizes L = 300, Lz = 500. 



5.2 Results 

In order to compare the analytical arguments presented 
in Section ^ and in the end to compute the total suscep- 
tibility of the manifolds we will first study numerically 
in (1 -|- 1) dimensions the finite size scaling of the aver- 
age lowest energy level. Then the shape of the energy gap 
distribution as well as the finite size scaling of the aver- 
age gap energies are considered. After that the numerical 
results of the first jump field distribution are presented to- 
gether with the finite size scaling of the average first jump 
field. The shape of the first jump field distribution and 
the finite size scaling of its average are shortly reported 
for (2 -I- 1) dimensional interfaces, too. The susceptibility 
is derived from the distribution and the finite size scaling 
of the first jump field. In the end of the section the jump 
distance distributions are considered. 

In order to see the logarithmic correction of the lowest 
energy level {E^), Eq. (^), when the height of a system, 
and thus the number of the minima, are increased, we have 
plotted in Fig. ||for three different lengths of the directed 
polymers L = 100, 200, and 300 



M-<^^[ln(^^.)]^/^ 



(25) 



where Nz ^ Lz/L'^, and {E) ~ 0.365L-I- 1. The prefactor, 
0.365, of the average energy of a polymer with only a single 
valley in a system, (i?), we have estimated by calculating 
systems of sizes Lz x L, where Lz = 6.5 x w is the aver- 
age roughness of a polymer, up to system sizes 50 x 1000 
with 2000 realizations. One should note, that the prefac- 
tor is highly sensitive to the disorder type and boundary 
conditions, c.f. Ref. and the constant factor (unity) 
also to the estimate of the size of the single valley. We have 



used it as a free parameter. Although for the following re- 
sults the relevant part of the tail of {E) should not be a 
pure Gaussian (at least if we have enough many minima 
Nzj so that the bulk of the distribution is avoided), we 
have used for 1/ry = 1/2 instead of l/r? = 1/1.6 = 0.625 in 
the fit of the logarithmic correction. In practice one can 
not tell these two choices apart in the range of the system 
sizes used. 

The probability distributions for the energy gaps of di- 
rected polymers in system sizes L^ = 100^ and 200^ are 
shown in Fig. ^(a). The distribution has a finite value at 
AEi — and the tail has approximately a stretched expo- 
nential behavior, exp(— aa;^), with an exponent h ~ 1.3. In 
the figure there is plotted as a comparison an exponential 
exp(— x) line, from which the deviation of the probability 
distribution is more clearly seen in the inset, where the 
distribution is in a natural-log scale. The finite value at 
AEi = is consistent with the weak replica symmetry 
breaking picture. 

To derive the scaling function for the AEi it is ex- 
pected that in systems with height Lz small enough to 
restrict the number of minima AEi mainly depends on 
the height of the system Lz- On the other hand, when 
Lz is large enough, there are enough valleys of which to 
choose the two minima, and one has AEi ^ AE ^ L^, 
hence 

(/i£;i(i:,i.))^|{l^^^' li'^l'^ (26) 

when the smaller parameter being varied. Since it is as- 
sumed, that Lz ~ i^, a natural scaling form based on 
these limiting behaviors is. 



(27) 



(AE,{L,Lz))-L''f[^ 



(28) 
I) by 



The argument y — Lz/L'^ for the scaling function f{y) is 
just a function of the number of the minima, i.e., Lz/ L^ ~ 
Nz^ and the scaling function has the form from Eq. (p4[), 
when ?7 = 2, 

f{y) ~ [lny]-i/2. 

In Fig. |^(b) we have plotted the scaling function (|2S_ 
collapsing {AEi{L, Lz)) / L^ versus Lz/L'' for various L 
and Lz, and find a nice agreement confirming the scal- 
ing behavior Eqs. (^^, ( pq ) as well as the analytic form 
Eq. ( p^ again assuming a Gaussian distribution (rj = 2). 

Next we explore the first jump fields. Consider the re- 
lation of the gap distribution and that of the jump fields 
given by Eq. (|l^). If we approximate P{AEi), Eq. (p^, 
with a uniform distribution, we get for the probability dis- 
tribution of the first jump field 



P{hi) — exp 



Nzhi 



hi 



In 



1 _ 111. 



1-hi 



exp(-/ii). 



(29) 

The form of P{hi) implies that it has a finite value at 
hi = 0, which is again consistent with the weakly bro- 
ken replica symmetry picture [p5| and also with the dis- 
cussion in Section M. We have also tried exponential and 
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power-law type of probability distributions for P{AEi) in 
Eq. (^3|), but all trials with negative exponents vanishes 
too fast with hi compared to the numerical data, and all 
positive exponents have behaviors with P(hi ^ 0) ^ 
and diverge for larger hi . 

In Fig. ||(a) we have plotted the probability distribu- 
tion of the first jump field for the system sizes = 100^ 
and 200'^. The probability distribution of the first jump 
field is similar to the probability distribution of the gap 
energies. The analytic formula, Eq. (p9|), is drawn as a line 
in the figure. One clearly sees that the line is a pure expo- 
nential, exp(— x), and the deviation of the numerical data 
from the exponential is similar to Fig. ^(a) of the energy 
gap distributions. Hence the shape of the numerical first 
jump field distribution is approximately a stretched expo- 
nential, exp(— ax''), with an exponent b ~ 1.3. Based on 
Figs.|(a) and|(a) one sees, that the probability distribu- 
tions of the energy gaps and first jump fields are the same. 
Thus with the correct P{AE) one gets from Eq. ([l^) the 
same P(hi) — P{AE). This distribution may be for small 
hi flat, but obviously starts to vanish for larger hi since P 
decays faster than exponentially. Another reason can be 
the fact that by discarding samples with the GS outside 
of a predefined window we just constrain the number of 
valleys in the sample so that the expectation value is the 
same in each one. can however fluctuate from sample 
to sample. The most important consequence is in any case 
that P{AEi = 0) = P{hi = 0) ^ 0. The finite size scaling 
and the normalization of the probability distribution of 
AEi give P{AEi = 0) - L-^[ln(L^L-^)]i/''. 

In order to find the scaling relation for the first jump 
field hi, we make the Ansatz (AEi) — {hi)LLz, since 
the field contributes to the manifold energy proportional 
to {D = 1) and Lz ^ i^zi) is the difference in the 
field contributions hz to the energy at finite h at different 
average valley heights zq, z\. Hence 



(30) 



where the scaling function /(y) for the number of the 
minima Nz ~ Lz/L'^ ~ y has the same scaling function 
Eq. (^). Fig. 11(b) shows the scaling function ( p8|) with a 
collapse of {hi{L, Lz))L^~^ Lz versus Lz/L'^ for various L 
and Lz with a good agreement, again. 

Next we move over to the (2-1-1) dimensional mani- 
folds. The inset of Fig. ^ shows the tail of the distribution 
for the system size L^ — 50^ with dilution type of disorder 
and bond occupation probability p = 0.5. The first non- 
overlapping jumps are included in the distribution, since 
due to the anomaly of the dilution disorder (lots of small 
scale degeneracy), there are typically two adjustments in 
the interface before the large jump. As a comparison, we 
plot again the exp(— a;) line, too. One sees that the devi- 
ation of the tail of the distribution from the exponential 
behavior is similar to the (1 -I- 1) dimensional case. The 
finite scaling of the first jump field in Fig. |^ is shown for 
interfaces of size L^ x Ly ~ 50^ and systems of height 
Lz = 30 — 90. We have fitted for a constant L the formula 



jOp, i.e., {hi{Lz)) ^ L^i[ln(Lz)]-i/2 and it works within 
error bars. 

Generalizing the numerical results of (1 -I- 1) and (2 -|- 
1) dimensional calculations and the analytic arguments 
from the previous section to arbitrary dimensions gives the 
behavior of {hi{L,Lz)) ~ L^-^Lji[ln(L JL9]-i/2. Since 
the probability distribution has a finite value at P{hi = 0) 
and (hi{L, Lz)) vanishes with increasing system size, one 
obtains from the normalization factor at P{hi = 0) for 
the scaling of the susceptibility, Eq. {m 



X ■ 



-D-e 



LzMLz/L'^)]'/^ 



(31) 



and in the isotropic limit, L oc Lz, the total susceptibility 
Xtot = L'^x becomes 



Xtot-L'^'+'-'ai-OHL)]'^'- 



(32) 



Thus the extreme statistics of energy landscapes shows 
up in the susceptibility of random manifolds or domain 
walls in a form that Eq. (|3^) has a logarithmic multiplier 
and the scaling behavior of the first jump field is due to 
the scaling of the energy minima differences. This is in 
contrast to Shapir's result and the logarithmic con- 
tribution is also missing from our earlier paper p3| where 
the valley energy scale was taken to follow the standard 
L^ assumption. Therefore one can conclude that the effect 
of extreme statistics is important in this problem: since we 
look at the finer details of the landscape "typical" differ- 
ences are not sufficient. Notice that for two dimensional 
random field Ising magnets C ~ 1 at large scales |]l^ and 
the susceptibility does not diverge |Q: the premise that 
A^z > 1 is broken in that case. If the condition A^^ > 1 
is violated the extreme statistics correction disappears. 
For flat interfaces so that the effective roughness expo- 
nent is zero (Ce// = 0), e.g. {100} oriented RB inter- 
faces with weak disorder and small system sizes L < Lc, 
when C = ||,|7[|38[||, 9 = D/2 and Eq. (H) becomes 
Xtot - i2D+i-D72[i^^(^)]i/2 _ L4[ln(L)]i/2^ when D^2. 

Finally we report the jump distance distributions of 
the first jumps. Fig. |lO|(a) shows P{Azi/zo) of the first 
jump with the field in (1-1- 1) dimensions. The distribution 
is a superposition of two behaviors: the interface jumps to 
the lowest minimum as such, and the external field favors 
the minima close to the wall (remember the differences 
of zi and zi* in Fig. Since the non-overlapping cases 
are excluded and the wall has a repulsive effect, there are 
cut-offs in the both ends of the distributions. The shape 
of the distribution does not change with the system size, 
which is consistent with the assumptions in Eqs. (^ and 
(|0l), {Azi) - Lz. Fig. [lO|(b) shows the same distribution 
of P{Azi/'zq) with a field for (2-1-1) dimensional manifolds 
and it is clear that the shape of the distribution does not 
depend on the dimension. In Fig. |lO|(c) we have plotted 
P{Azi/^o) without a field, i.e., the distance of the true 
lowest energy minima in (1 + 1) dimensions. Now there is 
no field, and thus the shape of the distribution is just a 
uniform one, again consistent with the predictions made 
in Section I in Eq. (|l|). 
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6 Scaling of the first jumps: field and distance valley to the waU we have 



Following similar arguments as Eq. in Section ^ for 
an excitation, a mean field result for the finite size scal- 
ing of the first jump field is next derived. Let us have 
an interface at an arbitrary height zq with an energy Eq 
and an isotropic system L (x Lz- The energy gap be- 
tween the two lowest energy minima, Eq. (M), scales as 
AEi - L^[{1 - C) ln(L)]~i/2, when the zq L^. However, 
when the interface is at an arbitrary height the number 
of the available minima depends on the original position 
of the interface, and we use only the dominating algebraic 
behavior so that AEi ~ . 

On the other hand the energy difference of elastic man- 
ifolds at different heights due to the field contribution is 
AE ~ hAzL^. Assuming that (Az) ^ L as in the earlier 
arguments, the field contribution becomes AE ^ hL'^. We 
expect that the first jump happens, when the gap equals 
the field energy, and thus the first jump field has a scaling 



{hi 



-e-d 



(33) 



In Figs. |lT|(a) and |ri|(b) we have plotted the aver- 
age first jump field {hi) in isotropic systems for (1 -I- 1) 
and (2 -I- 1) dimensional manifolds, respectively. The data 
is taken from non-overlapping jumps, to minimize finite 
size effects since they are smaller than with the overlap- 
ping jumps included. However, the fraction of overlapping 
jumps is small as seen in Fig. ^ for (1 -t- 1) dimensional 
case. The data in Figs. |ll| confirms within error bars the 
expected behavior of (hi) L~^/^ and {hi) L~^-^^ 
for (1-1-1) and (2 -f 1) dimensional manifolds, respec- 
tively. The logarithmic term [ln(A^z)]~^/^ might contribute 
a downward trend since it decreases with the number of 
the available minima N^, here disorder averaged due to 
the arbitrary GS location. This is however not noticeable 
in the data: the finite size effects in the calculated data are 
towards smaller absolute value of the exponent. In the in- 
sets of Figs. 0(a) and (b) the linear scaling of jump sizes 
{Azij ^ L is confirmed, which was assumed in Eqs. (^ 
and 1^). 

The evolution of the jump size distribution for succeed- 
ing jumps is demonstrated in Fig. |l^. As a comparison for 
Fig. ^ where only the non-overlapping jumps were con- 
sidered, here we have shown the overlapping ones, too. It 
is seen as a peak near Ziz„/z„_i = 0. The first jumps 
have clearly the most weight in the large jump end of the 
distribution, but for the following jumps the whole distri- 
bution shifts towards Azn/'Zn-i = 0. After a small number 
of jumps the interface is already in the minimum closest 
to the wall and then the random-bulk wetting behavior 
takes over. In order to estimate the number of the jumps 
the manifold does before binding to the wall, it is assumed 
that after a jump the system looks statistically the same 
as before it. From Fig. [lO|(a) we infer that the probability 
distribution for the jump size has a form P{Az) ~ Az'^ 
(note that we did not attempt to compute this from the 
analytical valley arguments) and that is then taken to be 
the same for all (large) jumps Azn. To calculate the aver- 
age jump size of jumps which do not jump to the closest 



{AZn) 



Zn-l—Ai W 



Bn^iAziAz)^ d{Az) 
_ 3(Zn-i - Aiw)'^ 



(34) 



where the upper bound of the integral is to neglect the 
cases that jump to the closest valley to the wall, Ai is the 
prefactor to multiply the roughness value to get the val- 
ley width, and i?„-i ~ "ij z\_i normalizes the probability 
distribution. Using (z„) = (z„_i) — {Azn) we get the next 
position of the interface. In order to calculate the proba- 
bility of an interface to jump with the first jump to the 
closest valley to the wall, pi, we use the same probability 
distribution as in Eq. 



Pi = / BoiAz)^d{Az) = 1- 



1 - 



Aiw 



= Pi, 



and similarly for the probability of an interface to jump 
with the second jump to the closest valley to the wall 



P2 ^ {I ~ Pi)p2 = [l ~ Pi) Bi{AzYd{Az) 

'J zi —Aiw 



1 - 



Aiw 



1 - 



1 - 



Aiw' 
zi 



Pn = Pr 



where 



fe=l 



-Pk) 



(36) 



l-{l-qnf\ n(l-*)'' (37) 



Due to the hierarchy we finally get for the n jump 



fe=i 



Aiw 

Zn-l ' 



(38) 



To get an estimate for the number of jumps, we write 



npn- 



(39) 



This can be evaluated, but only approximately since among 
others the estimate of Eq. (^^ breaks down beyond n 
smah. Using Aiw ~ AL'^ = 50, L = 1000, which is the 
case with Ai ~ 6.5, (see the numerics in Section 5.2), 
and taking zq = 1000, we get (A/") ~ 3. Eq. (||) gives 
for the first jump size, with the above numerical values, 
{Azi) ~ 0.620, which is not far from the behavior in the 
inset of Fig. |ll|(a), where {Azi) ~ OAL^ ^ 0.8zo, smce 
(zo) - Lz/2. 



7 Application of a non-zero field: 
random-bulk wetting 

In order to see the wetting behavior one just studies the 
effect of a large external field on the average interface to 



12 



E. T. Seppala and M. J. Alava: Energy landscapes, lowest gaps, and susceptibility of elastic manifolds at T = 



wall distance. Notice that the low-field physics discussed 
extensively above can be considered as the eventual out- 
come in any sample with Lz > L'^ so that more than one 
valley is available. This means that when the field is de- 
creased from a large field value, the interface finally jumps 
into the bulk. Fig. |l^(a) we show the average mean height 
{'z{h)) versus the field h for (1-1-1) and (2-1-1) dimensional 
manifolds. The system sizes are x L — 100 x 3000 and 
LzXL^ = 50 X 300^. To maximize the prefactor, A2, in the 
scaling of roughness, w ^ and hence the width of 

the minimum energy valley, the disorder has been chosen 
to be strong, i.e., dilution type of disorder with small p, 
but above the bond-percolation threshold. However, there 
are still some deviations in the form of greater exponents 
than the expected from Eq. (|^), which gives the values 
ip = 1/2 and ip ~ 0.26 in (1 -I- 1) and (2 -I- 1) dimensions 
from C(i-i-i) — 2/3 and C(2-i^i) = 0.41 ± 0.01, respectively. 
We found that the trend is nevertheless clear, with greater 
L and fixed Lz the exponents become closer to the ex- 
pected one. The effective exponent ipeffiL) can be used 
to extract the asymptotic, L-independent exponent in par- 
ticular in (2 -|- 1) dimensions since this case is hampered 
most by finite-size effects. ipeff{L) as a function of 1/L 
indicates that the asymptotic value is indeed 0.26 ± 0.02 
and that the system sizes at which ipeff{L) approaches 
that are of the order of L = 10"*. That is, only with 10^" 
sites it becomes possible to reach the asymptotic regime. 
When one has L oc Lz calculating the average mean height 
{z{h)) with a fixed field h is nothing but averaging over 
the jumped and not jumped interfaces together with their 
location, see Fig. ||(a). 

In Fig. |l3|(b) we have plotted as a comparison the av- 
erage mean height {'z{h)) versus the field h for (2 -I- 1) 
dimensional manifolds with a weak disorder. In this case 
the weak means for system sizes used, that the rough- 
ness of the manifolds are not in the asymptotic rough- 
ness limit yet, L < L^ |^,|^,|8[^. The system sizes 
used are x = 50 x 80^^— 400^ and the behavior 
is simple: either the interface stays in its original posi- 
tion or jumps directly to the wall. Taking into account 
the jumped and original interfaces as {z{h)) = {z[l — 

P(zo, /i)]^o + P(5o, h)xO) = {[l- P{zo, h)]zo) = /o''-- [1 - 
P{'zo{h)}]'zo dzQ, gives the behavior {z{h)) ~ i.e., the 
exponent ^ ~ 1, if P{zo(/i)} ^ h~^. The larger mani- 
folds jump faster to the wall, i.e., they feel the pertur- 
bation earlier, since Eq. (pG) for flat interfaces (C = 0) 
becomes {hi{L,Lz)) - L''^Lji[ln(L^)]-i/2. with fixed 
Lz and 6 — D/2 from Poissonian statistics, we get in 
D — 2 {hi{L)) L^^ consistent with the numerical 
data. This leads to the behavior of the wetting scaling, 
{z{h)) ~ c{L)h-'f' where c{L) ~ L"! and V = 1- The 
finite size scaling of the prefactor indicates that at large 
I/-limit with fixed Lz the flat interfaces are immediately at 
the wall, and thus the systems are non-wet. This implies 
that there is an interesting cross-over around Lc between 
such a "dry" regime and the bulk wetting that takes over 
for still larger L.ln D > 2 this kind of behavior is relevant 
even in the thermodynamic limit, if the disorder is weak 



due to the presence of a bulk roughening transition for 
bond disorder. 



8 Discussion 

8.1 Finite temperature behavior 

The movement of the elastic manifolds in random media 
at low temperatures, when an applied force is much below 
the depinning threshold Fc, is characterized as creep. The 
dynamics is controlled by thermally activated jumps over 
pinning energy barriers, which separate the metastable 
states. D. Fisher and Huse ||l^ showed that for a DP at 
finite temperature T > the fluctuations of the entropy 
(Z\S')^ and the internal energy (AEint)'^ scale linearly with 
the length of the polymer and cancel each other. Hence 
there are only the fluctuations of the free energy (Z\F)^, 
which scale with the zero temperature energy fluctuation 
exponent 26 — 2/3. Since the free energy is the one which 
should be minimized at flnite temperature, it is the one 
which defines statistically the shape of the energy land- 
scape, although the energy valleys and minima need not 
to have exactly the same real space structure as at T = 0. 
Thus the 6* = 1/3 exponent should define the energy gaps 
also when T > 0, expect in the cases there is a critical 
temperature Tc- Hence, our derivation of the susceptibil- 
ity and also the first order character in the reorganization 
of valleys should be relevant also at T > 0. 



8.2 First arrival times in nonlinear surface growth 

The (l+n) dimensional directed polymers map, in the con- 
timmm limit, to the KPZ ||,^ equation by associating 
the minimum energy of a DP-configuration with the min- 
imum arrival time ti = i?o of a KPZ-surface to height H. 
The connection is illustrated in Fig. |lj. The minimal path 
of the DP with the end point Xi(ii) equals the path by 
which the interface reaches H, at location xi and at time 
ti. Thus ti attains a logarithmic correction, from Eq. (p2[), 
of size -H^{\n{L/H^/^)Y/'i, where L is the linear size of 
the system, P = 9 > and z = 1/C are now the roughening 
exponent and dynamical exponent of the KPZ universal- 
ity class and the values of 6* = 2C — 1 and ^ > 1/2 depend 
on n. Notice that if there is a upper critical dimension Uc 
in the KPZ growth, then the logarithmic correction is not 
there anymore, and 9 = 0, = 1/2, i.e., a random walk 
ensues. Consider now the second smallest arrival time ^2- 
In the language of directed polymers, if the path X2(i') 
corresponding to t2 is independent of Xi(i') that corre- 
sponds to ti the time and the path are found inside a sep- 
arate, independent valley. The difference At — t2 — ti then 
equals AEi, and likewise obeys extreme statistics, so that 
At - i/^[ln(L/i?i/^)]-(''-i)/''. For growing surfaces this 
limit is the early stages of growth, in which the correlation 
length ^ ^ L, and therefore the arrival times, or directed 
polymer energies, are independent. On the other hand if 
we disturb the growing process such a way that it depends 
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on the X, e.g. linearly with a factor h |Q , the polymer end- 
ing at X2 becomes faster ii h > hi, see right hand side of 
Fig. Similarly now the factor hi has a scaling behavior 
from Eq. @), hi - i/^'-iL-i [ln(i/i?i/^)]-(''-i)/''. 



9 Conclusions 

To conclude we have studied the (1 + 1) and (2 + 1) dimen- 
sional elastic manifolds at zero temperature, when an ex- 
ternal field is applied. We have demonstrated that the re- 
sponse of manifolds shows a first order character ("jump") 
in the sense that the manifolds change their configuration 
in the large system size limit completely. This persists in 
a finite system over a small number of such jumps. The 
distance that the center-of-mass moves is extensive. The 
whole picture is based on a level-crossing between two low- 
energy valleys in the energy landscape. Averaging over 
the total magnetization of such random magnets with a 
domain wall or over the positions of manifolds when an 
external field is applied becomes dependent on whether 
the DW or manifold has jumped or not. This leads to 
the problem of self-averaging in random systems. Here a 
disorder average smooths over the "coexistence" between 
systems that are not affected by a finite field h and those 
that have responded. 

In order to study the susceptibility of the DW in ran- 
dom media, one has to take into account the probability 
distribution of the sample-dependent field associated with 
the global change of the configuration, and take the limit 
of vanishing fields. This probability distribution has a fi- 
nite density at h = 0. The finite size scaling of the density 
is dependent on the finite size scaling of the number of the 
low lying nearly degenerate energy minima in the system. 
The scaling of the number of the energy minima leads to 
a logarithmic factor in the susceptibility, and can be ac- 
counted for by using extreme statistics. Such effects are 
difficult to study by usual field theoretical means, since 
one has to have access to the whole probability distribu- 
tion and not only the few first moments thereof. Notice 
that the crucial difference to much previous work is the 
simple fact that we allow for multiple minima in the energy 
landscape, which is often most excluded by the boundary 
conditions applied to the problem. Although the deriva- 
tions and the numerical calculations done here have con- 
centrated on random bond type of randomness similar be- 
havior should be seen in random field cases, too. 

The discrete character in the movement of elastic man- 
ifolds with an external field results also in that the con- 
tinuum theory for wetting in random systems works only 
in slab geometries, where there is room only for a sin- 
gle valley, or in the large external field limit, when the 
interface is close to the wall. On the other hand the fiat 
interfaces are shown to jump directly to the wall, i.e., to be 
non-wet. It would be interesting to see if the dynamics of 
the manifolds at finite temperature reflects the flrst order 
character seen here at T = 0. Through the connection of 
(1+1) dimensional DW's to the directed polymers of the 
KPZ surface growth, we have shown that the logarithmic 



factor is also present in the statistics of growth times in 
nonlinear surface growth. 
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Wj 



W2-Z2 



Fig. 1. A schematic figure of two interfaces in a system of 
parallel length L and perpendicular one Lz > L'' at heights zi 
and 22 — W2, with the corresponding roughness values wi and 
'W2- C-L denotes the transverse correlation length of a part of 
the interface of length . See the text. A droplet is seen in the 
interface of height zi as a dotted line. 
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Fig. 2. (a) Examples of two realizations of changes in mean 
heights 2 of interfaces normalized by their original (global min- 
imum) positions zo vs. applied field h for (1 + 1) dimensional 
systems. The change in the external field is done in steps of 



Ah = 10 . Note the large jumps in both cases. L 



200" 



yEi(h)=Ei+hL"zi 

E2(h)=E2+hL% 




Fig. 3. The level-crossing phenomenon for interfaces in ran- 
dom systems in the presence of an external field. Originally the 
interface lies at height zo and has an energy Ea . When the field 
is applied its energy increases linearly and at h-i the interfaces 
jumps to 2i < 00, with the energy Ei{hi) = Ei + h\L^ zi, 
where Ei is the energy of the interface at zi without the field. 
Similar behavior takes place at /),„, n = 2, 3, . . ., when the inter- 
face moves from Zn-i to z„ < Zn-\- The thick line represents 
the minimum energy of the interface E(h) as a function of the 
field h, with discontinuities of the derivative at 

„ -1.40 

3 



> 
o 



200 .-ti 



C/3 



-1.60 



-1.80 



-2.00 



-2.20 




Jij.z £ [0 — 1] uniform distribution and Jij^x = 0.5 (random 
bond disorder) . (b) The expected scenarios (droplet formation, 
jump to the lower edge of the system) before and after the first 
moves from global minima zo{x) to zi{x). 



3.5 4.0 4.5 5.0 5.5 6.0 6.5 

system size [ln(L)] 



7.0 



Fig. 4. The fraction of jumps of interfaces with initial global 

minimum in an arbitrary position that have an overlap between 
interfaces before and after jump q as function of the system 
sizes L. The system sizes are = 50^ — lOOO'^ and the number 
of realizations N = 3000 — 5000. The line is a least-squares fit 
to the data with the power-law behavior 7^-o-23±o.oi 
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Zo z 




^ Mj,^- kzo-hZj* 
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Fig. 5. (a) A simplified view of the minima in the energy 
landscape of a random system. AEi = Ei — Eq is the energy 
gap between the ground state at zo, denoted with a black circle, 
and the second lowest minimum at z\ denoted with a gray 
circle. AEi* = Ei* — Eq is the energy gap (normalized with 
L^) between the ground state and the minimum at zi* . (b) 
The view of the minima in the random system when the field 
is applied. At hi the interface moves from to z-i* , indicated 
with a gray circle, since the energy difference AEi* = Ei* — 
Eq = hizo — h\z\* = while all the other AE^s are greater. 
However, often E\' and E\ coincide. 
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transverse system size (L^/L ) 

Fig. 6. The scaling of the ground state energy, ((_E) — {E(,))/L^ 
vs. Lz/L''', see Eq. ( |2^ ) and the text, for the system sizes 
L =100, 200, and 300, each with zo/L^ ~ const. 0=1/3, 
C = 2/3. The line -0.41 + 0.53[ln(2.78L3/L« is a guide to 
the eye. The number of realizations ranges from A'^ = 500 for 
L = 300, Lz = 500 to A = 2000 for L = 200, = 600. 
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Fig. 7. (a) The scaling function of the probability distribution 
P{AEi/{AEi}) X (AEi) for the energy gaps AEi of the lowest 
two minima normalized by their disorder-average {AEi) in a 
(10-base) semilog-scale for the system sizes LxLz = = 100^ 
and 200^ . The inset shows the tails in the natural-log-scale. The 
initial global minimum position zq/L^ ~ const for all L. The 
number of realizations N — 2 x 10* for both system sizes. The 
lines in the figure and in the inset are exp(— i). (b) The scaling 
function f{y) of the scaled disorder-average of the energy gap 
{AEi)/L^ as a function of scaled transverse system size L^/L^ 
for the system sizes L =100, 200, and 300, each with zo/L^ — 
const, e = 1/3, C = 2/3. The line f{y) = 0.23ln(y)-^^^ is 
a guide to the eye. The number of realizations ranges from 
= 500 for L = 300, = 500 to iV = 2000 for L = 200, 
Lz = 600. Only the non-overlapping interfaces are considered. 
The fraction of overlapping ones is vanishingly small. 




I 



-3.0 -2.0 -1.0 0.0 1.0 

normalized jump field lg[hj/<hj>] 

0.50 p 
0.45 



^ 0.40 
V 0.35 



0.30 



0.25 
.3. 0.20 




o L=100 
□ L=150 
o L=200 
A L=250 
< L=300 



,U 0.15 



10" 



10' 



transverse system size (L/L ) 

Fig. 8. (a) The scaling function of the probability distribution 
P{hi/ (hi)) X (hi) for the first jump field values hi normalized 
by their disorder-average (^i) in a (10-base) semilog-scale for 
the system sizes Lx Lz = L"^ — 100^ and 200^ . The inset shows 
the tails in the natural-log-scale. The initial global minimum 
position zo/Lz ~ const for all L. The number of realizations 
N — 2x 10* for both system sizes. The line is the analytic result 
from Eq. (|l^) with a uniform distribution P{x) and Nz = 20. 
(b) The scaling function f{y) of the scaled disorder- average of 
the jump field {h\)L^~^ as a function of scaled transverse 
system size LzjL'^ for the system sizes L =100, 150, 200, 250 
and 300, each with zo/Lz ~ const. 9 = 1/3, C, = 2/3. The 
line f{y) = 0.41 ln(t/)~^/^ is a guide to the eye. The number 
of realizations ranges from TV = 500 for L = 300, = 500 to 
N = 2600 for L = 200, = 600. Only the non-overlapping 
jumps are considered, again. 
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Fig. 9. The disorder-average of the first jump field (hi) as a 
function of transverse system size Lz for the (2+1) dimensional 
interface. The disorder is random bond type, with SJij/Jo = 1, 
the number of realizations for each system size A'' > 100, and 



the system size Lx 



50. The initial global minimum po- 



sition zo/Lz ~ const. The line 0.015I/J^[ln(I/z)]~^/^ is a guide 
to the eye. The inset shows the tail of the scaling function 
of the probability distribution P(h\/ {hi)) x {h\) for the first 
jump field values h\ normalized by the disorder-average (/ii) 
in a natural- log-scale for the system size = 50^. The initial 
global minimum position zq/Lz — const. The disorder is dilu- 
tion type with p = 0.5 and the first non-overlapping jump is 
considered. The number of realizations N > 3000. The line is 
the same as in the inset of Fig. ra(a). 
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Fig. 10. (a) The histogram P{Azi/zo) of the jump size Azi 
of the first jump normalized with the initial position zq (the 
initial global minimum position zo — const) for (1 + 1) dimen- 
sional interface. The number of realizations N = 2 x W^, the 
system size — 100^, and the random configurations are the 
same as in Fig. ^a). (b) P{Azi/zo) for (2 -I- 1) dimensional 
interface The number of realizations N ~ 3000 and the ran- 
dom configurations are the same as for the data in Fig. |^ inset. 
The system size = 50^ and disorder is dilution type with 
p = 0.5. The initial global minimum position zo — const, (c) 
P{Azi/ Zq) for (1 -|- 1) dimensional interface without a field, 
= 100^, the number of realizations N = 2 x lO**. The ran- 
dom configurations are the same as in Fig. |^a) and the initial 
global minimum position zo — const. Note in all figures the re- 
pulsion induced behavior at the both ends: for values close to 
zero due to neglect of the overlapping interfaces and for values 
close to unity due to wall repulsion. 



E. T. Seppala and M. J. Alava: Energy landscapes, lowest gaps, and susceptibility of elastic manifolds at T = 19 



10" 



10 



^ _ 

•:?io 

t/3 

10" 



10" 



A 
V 



y:^ 10 



I 

J3 



C/5 



10" 



1.62±0.04 




200 400 600 800 1000 
system size (L) 



10^ 



10^ 



10 



system size (L) 



(b) 



6.56L 



-2.18 




50 60 70 80 
system size (L) 



10^ 



10^ 



system size (L) 



Fig. 11. (a) Finite size scaling of the average first jump field 
(hi) for DP's. The line is the least squares fit to data. The scal- 
ing argument gives a = 5/3. The inset shows the average jump 
distance at the corresponding field hi with a linear fit 

to data. ( ) is the disorder-average over N = 1000 realizations 
for the system sizes L x = L'^ = 50^ and 100^, A'' = 500 
for = 200^ - 400^ and = 200 for = 600^ - 1000^ 
The disorder is of random bond typo. No overlaps included, 
(b) Finite size scaling of the average first jump field (hi) for 
(2 + 1) dimensional manifold. The line is a guide to the eye 
with a slope a ~ 2.18. The inset shows the average jump dis- 
tance (Azi) at the corresponding field hi with a linear fit to 
data. ( ) is the disorder-average over N ~ 250 realizations for 
each system size. The disorder is dilution type with p = 0.5 
and the first non-overlapping jump is considered (usually two 
small-scale adjustments take place before the jump). 
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Fig. 12. Evolution of the histogram of the jump sizes 
Az„/z„-i normalized with its previous position for n — 
1,2,3,4. The system size is = 100"^, the number of real- 
izations N = 2000. The overlapping jumps are included. 
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Fig. 13. (a) The average interface mean height {z{h)), i.e., the 
thickness of wetting layer, as a function of the external field h 
for one dimensional directed polymer, open circles. The system 
size is L = 3000, Lz = 100 and the number of realizations N = 
1000. The disorder is dilution type with p = 0.55. The filled 
diamonds denote (2 + 1) dimensional interfaces. The system 
size is Lx ~ Ly — 300, Lz = 50 and the number of realizations 
A'^ = 30. The disorder is dilution type with p — 0.30. The solid 
line is a guide to the eye with a slope ^p = 1/2 and the dashed 
line has a slope ip = 0.26. (b) (z{h)) vs. h for various system 
sizes Lx = Ly = L = 80, 100, 120, 200, 300, and 400 in the flat 
regime. The disorder is dilution type with p = 0.95. Lz — 50 for 
each system. The number of random configurations is Af = 100- 
- 200 for L = 80 - 120, and AT = 30 for L = 200 - 400. The 
line is a guide to the eye with a slope ip = 1. 
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Fig. 14. The relation between directed polymers and growing 
interfaces. Two directed polymers in independent valleys equal 
the fastest arrival time ti at xi, solid line, and the second 
fastest at X2 with time 1-2, dotted line, of a KPZ interface to 
a fixed height H , when the external field h = 0. In the right 
hand side figure an external field is added, which increases the 
growth time depending on the position in the direction of the 
arrow and the polymer at X2 becomes the one corresponding 
to the fastest time to reach H. 



